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Linear guided waves in hyperbolic slab waveguide. Dispersion relations 
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Guided waves modes in a slab waveguide formed from the isotropic dielectric layer embedded by 
hyperbolic materials are investigated. Optical axis is normal to the slab plane. The dispersion 
relations for TE and TM waves are found. The differences between hyperbolic waveguide and 
conventional one are demonstrated. In particular, for each TM mode of hyperbolic waveguide there 
are two cut-off frequencies and the number of modes is limited. For the TE and TM modes Poynting 
vector component along the wave’s propagation axis could be equal to zero. 
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I. INTRODUCTION 

Metamaterials and their optical properties attract grate attention over the past decade. Metamaterials are artifi¬ 
cial materials which are composed of unit cells far below the size of the wavelength. These materials can exhibit 
exotic electromagnetic properties. The negative refraction is the famous property of these media CHS]. As usually 
such materials are materials with simultaneously negative real parts of the dielectric permittivity and the magnetic 
permeability in some frequency region lain. 

It has been known, that the negative refraction can take place in the anisotropic media HCD. It should be pointed, 
that uniaxial anisotropy is typical property of the metamaterials. It is important to remark, that an anisotropic meta¬ 
materials can demonstrate negative refraction in one direction and positive refractions in the orthogonal directions. 

Let us assume, that in the uniaxial anisotropic medium coordinate axes OX, OY and OZ are chosen to be equal to 
the principal axes of the dielectric permittivity tensor with the principal values of dielectric permittivity, which satisfy 
the conditions Sxx = and Syy = = Eq- The dispersion relation for extraordinary wave connecting frequency uj 

with Cartesian components of wave vector k, takes the following form: 

ee{uj) So{oj) c2 

This relation shows, that in the case of either or is negative, iso-frequency dispersion surface 0 represents the 
hyperboloid. The hyperboloid of one sheet is realized if Ce > 0, Sq < 0 and hyperboloid of two sheets is obtained if 
Se<0,So>0 [HHH]. 
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These anisotropic materials are referred to as hyperbolic materials. 

The hyperbolic materials can be fabricated as a multilayer structure consisting of alternating metallic and dielectric 
layers [muni [201, or as a nanowire structure consisting of metallic nanorods embedded in a dielectric host imi2iii22i- 

Iso-frequency dispersion surfaces 0 allow the infinitely large wave vectors. It results in the different effects, among 
which are the Purcell enhancement of the spontaneous emission rate in hyperbolic metamaterials [231126] and the 
subwavelength resolution effect [HEt]. 

The optical phenomena on interface between conventional dielectric and hyperbolic metamaterial has attracted 
attention. The surface waves were studied in [28]. The extremely large Goos-Hanchen shift has been studied in 
some details in [29]. Plasmonic planar waveguide cladded by hyperbolic metamaterials (Eig.l) was proposed and 
investigated in (201 Ell. 
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The purpose of this paper is to investigate the dispersion properties of the linear guided waves in a planar waveguide 
previously studied in [30]. Unlike [30] the guided waves localized in the dielectric core will be considered here instead 
of surface waves. The anisotropy axes of substrate and cladding layers are directed along OX axis that is normal 
to the interface (Fig.l). In the planar geometry, as is known [32], the Maxwell equations can be separated into two 
uncoupled systems of equations, which describe the propagation of the waves having a different polarization. These 
waves are refereed to as TE and TM waves. Analysis of the guided TE and TM waves will be done independently. In 
both cases the expression for dispersion relation connecting effective waveguide index with frequency will be obtained. 
The symmetrical waveguide will considered in details. 


II. ELECTRIC AND MAGNETIC FIELD DISTRIBUTIONS FOR GUIDED WAVES 

Let us consider a slab waveguide. We assume, that the material of the waveguide core is nonmagnetic /i^ = 1 and 
has an isotropic permittivity Si. The core thickness is h. The dielectric core is cladded by the uniaxial hyperbolic 
metamaterials which are characterized by symmetric dielectric tensor with the principal dielectric constants , 

So\se and the magnetic permeabilities jui and /is. All permeabilities are assuming to be positive. The anisotropy 
axes are aligned with a unit normal vector to the interface, i.e., along the OX direction. (Fig.l). Axes OY and OZ 
are parallel to interface. Axis OZ is directed along the wave propagation. In this case the Maxwell equations are 
invariant under the shifting along OY axis. Thus, the strengths of the electric and magnetic fields of the guided wave 
are independent of variable y. From it follows, that the Maxwell equations are splitting into two independent systems 
of equations describing the TE and TM waves [32]. 



EIG. 1: A schematic illustration of the hyperbolic waveguide 


The TE wave is defined by the tangent component of the electric field vector Ey and by two components of magnetic 
field: and Hz- These values are assumed to be harmonic functions of the time: exp(ico’t). The wave equation for 

complex amplitude E = Ey{x^ z^uj) looks like: 

^ ^ fcoeo(a:)M(a^)-E = 0, 

where ko = cj/c, uj is frequency of radiation. The principal dielectric constants and magnetic permeabilities are 
piecewise functions (Eig.l): 


So{x) = 


The components of magnetic field can be found from the following relations: 

i dE i dE 


to 

X < 0, 


X < 0, 


X < 0, 


0 < X < h. 

Seix) = < e* 

0 < X < h. 

k'ix) = < /ii 

0 < X < h. 


X > h. 


X > h. 

[ /^3 

X > 


= 


koy{x) dz ’ 


H. = - 


kQii{x) dx 


(3) 


The TE wave is ordinary one according to choosing of the direction of anisotropy axis. In the case of 5o > 0 the 
problem is reduced to a familiar case. However, one can expect the interesting result for hyperbolic material with 
< 0 . 
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The electric field can be represented in the form E{x^z) = exp(i/3z), because the waveguide is homogeneous 
along OZ, where parameter p is the propagation constant. This parameter is similar to the wave number in the case 
of homogeneous medium. Solutions of the wave equation should be found with taking into account the boundary 
conditions E^O, H^Oat \x\ oo. The solution can be obtained by standard procedure [32]. Distribution of the 
electric field is given by the following equations 

x<0: Ec.c., 

0 < X < h : E^‘^^ = A [cos(/^x) + sin(/^x)] + c.c., (4) 

X > h : E^^^ = A [cos{Kh) + ^p sin(/^h)] ^-Ax-h)^i/3z ^ ^ 

where following parameters are used 

+ fcoAi3|4^^|, = kluiEi - 0^. 

These parameters can be used to define the Goos-Hanchen phase shifts (f)q and <pp: 


^ = -tan((/)g/2) = ^p = -tan(0p/2) = 

K/lS K/J^i 

Normalized amplitude A of the electrical field at x = 0 is arbitrary. 

TM wave is determined by the component of a magnetic field Hy and an electric field components E^^ E^. The 
wave equation for magnetic field H = Hy{x^ z,cj) is 


1 d^E d 


1 dH 


+ kQ/j.{x)H = 0. 


Se{x) dz‘^ dx \So{x) dx 
The electric field components can be found from the following expressions: 

i dH ^ i dH 


Ex — — 


koSeix) dz ’ 


E,= 


koeo{x) dx ' 


(5) 

(6) 


The principal dielectric constants and permeability are piecewise functions previously defined in the case of TE waves. 

A consideration of the equation ^ shows, that if < 0 and > 0 (a = 1, 3), then there is no solution of this 
equation decreasing at |x| tends to infinity. Hence, there are no wave localized in waveguide. Alternatively, if > 0 
and S0'’ < 0 then under conditions; 




(7) 


equation admits a solution describing wave confinement in this waveguide. The magnetic field distribution is 
written as 


X < 0 : + c.c., 

0 < X < h : = A [cos(«;x) — sin(«;x)] + c.c., 

X > h : = A [cos(Kft.) — sin(Kft.)] + c.c. 


(8) 


In these expressions following parameters 


u(l) I , , u(3) I 

= klfii\e^^^\ -= kliiiSi- 


d3) 


are used. The GoosHanchen phase shifts (pq and (pp are defined by the following expressions: 


^q — tan(0g/2) — 



^p — tan(0p/2) — 


psi 
(i)i 
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III. DISPERSION RELATIONS 

Taking into account decreasing of the magnetic and electric fields in the limit \x\ oc, the Maxwell equations 
solutions describe the waves which are confined by the waveguide. It is necessary to distinguish the coupled surface 
waves and the waveguide modes. In a linear waveguide the amplitude of surface wave takes the maximum at interface. 
In the case of the dielectric waveguide cladded by metal this coupled surface wave is said to be the plasmon-polariton 
wave. Guided plasmon-polariton wave in the dielectric waveguide arranged from a hyperbolic metamaterial has been 
proposed and studied in [30l |3T] . The dispersion relation of this plasmon-polariton wave can be obtained by changing 
parameter to k? = — kliJ^iCi. However, other than these, there are a number of waves localized in the dielectric 
core that are designated as guided modes or waveguide modes [32]. 

For all guided waves the dispersion relation exists as expression connecting the propagation constant (3 and frequency 
uj. The dispersion relations of the guided waves under consideration can be determined by using the continuity 
conditions of an electric and magnetic field on the interferences. It is suitable to get the dispersion relations for TE 
and TM wave separately. 


Case of TE wave 


Distribution of the magnetic field in waveguide is derivable from 0 with taking into account the expressions 0. The 
continuity conditions for tangent components of both electric and magnetic field vectors result in following relation: 


^2iKh 




i + iij Vi + ^Ci 


= 1. 


Using the expression for the Goos-Hanchen phase shift one can write the dispersion relation in form: 

2i<ih + <pp + <pq = 27Tm^ m = 0, 1, 2,... (9) 

If the effective index of refraction rig/ is defined according to formula p = koUef^ than equation ^ can be written as: 


hkn 


nj — = arctan f ^ 




arctan f ^ .t 
V/^3 


/^3+^e/ 


(3), 


nf = jXiSi 


are introduced. 


Here the indexes of refraction nl = 

The dispersion relation shows that the effective index of refraction is limited by the condition 0 < < nj. In the 

case of all dielectric waveguide the similar limitation appears as max(ni,n 3 ) < < nj. Difference between these 

inequalities is due to the fact that in dielectric waveguide embedded in a hyperbolic media the total internal reflection 
takes place at any incident angle. As in the case of dielectric waveguide cladded by metal. 

Analysis of the dispersion relations will be performed for the case of symmetric waveguide, where nj = n\. Then 
the dispersion relation takes the form: 


koh 


nj — njj^ = 2 arctan 



nm. 


( 10 ) 


Equation (10) can be rewritten in the normalized form. The parameter b is introduced by the following equation 


nj-\-nj^ = 6A, where A = nj-\-nj. The normalized waveguide thickness V is introduced by formula V = k^h^Jnj 
Then the relation (10) takes the form: 


V \/l — b = 2arctan { 

\ /ii V l-b 


( 11 ) 


This equation defines the function 6(U, m) that is implicit dependence of the normalized effective index of refraction 
b on normalized waveguide thickness V. In the case under consideration b is in the interval [6o,l), where bo = 
njl{nj + nj). The plots of b{V,m) vs V are presented on Eig.2 at = 1-2 and bo = 0.2. The fact that 6o > 0 

means that zero-mode TEq has non zero cut-off frequency Ugo- Substitution b = bo into equation 0 results in 


UgO = 2W1 H—^ arctan 


/ lijni 


^cO — - A/ - I 2 \ 


( 12 ) 


Eor conventional dielectric waveguide cut-off frequency Vco is zero. 
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FIG. 2: Dispersion curves for TE modes of the hyperbolic waveguide. 


B. Case of TM wave 


By the use of (§ and (H the electric field strength cab be found. The continuity conditions for tangent components 
of both electric and magnetic field vectors lead to the following dispersion relation for TM waves 


^2iK,h 




V 1 - 


1 . 


Using the Goos-Hanchen phase shifts the dispersion relation can be written as : 

2nh-\-<pp(^q = 27vm^ m = 0, 1, 2,... 


In terms of initial variables the relation (13) takes the form: 


(13) 


h\/kQ{jj.iei — = — arctan , 




d3) 






(3)u(3) 


— arctan / 


- n 


ef 


el 


\ ko 


dl) 


^e/ 


(i)u(i) 


- n; 


ef 


Let us consider only symmetrical waveguide, where The dispersion relation can be rewritten as 


hko 



—2 arctan 


el A ^e/ 


+ 7rm, 


(14) 


where the effective indexes of refraction rii and Ue are used. These parameters are defined by correlation nf = fiiSi 
for isotropic dielectric and = /ii^e for extraordinary wave in the hyperbolic media. 

The condition for effective index 


0 < nlf < n,, 0 < nif < n\ 


- '^ef ^ 


follows from the equation ( |14[ ). In the case of a convenient dielectric medium this condition appears as ric < n^f < 
where ric is the refraction index of the substrate or cladding layer. 

Equation (14) can be rewritten in terms of the uniform variables (6, V) by using following correlations: = 

bA > 0, where A = n?—and normalized waveguide thickness V is introduced as U = kohy/nl — 71^. b is normalized 
effective index of refraction of waveguide. It results in the dispersion relation in following form: 



V-\-b = — 2 arctan 


+ irm. 


(15) 
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In the case under consideration the parameter b is in interval [0, 6o], 
The dispersion curves corresponding to equation (15) are represented 
6o = 2 is hold. For comparison the dispersion curves corresponding all 


where bo = (nf — nl). 

in Fig.3. The assumption = 1.2, 

dielectric waveguide are shown in Fig.4. 


V — b = 2 arctan 



7rm, m = 0,1, 2,... 


where u = ell{eoSe)’ The curves in Fig. 4 were obtained dX u = 1.2. 



FIG. 3: Dispersion curves for TM modes of the hyperbolic waveguide. 



FIG. 4: Dispersion curves for TM modes of the conventional dielectric waveguide. 

Figures show, that in the case of TM wave the number of guided modes of the hyperbolic waveguide is limited. As 
the dielectric core thickness h increases, one mode disappears but other mode appears. Furthermore, the zero-mode 
{m = 0) in this waveguide is absent. For an all dielectric waveguide (Fig.4) the number of guided modes increases 
with core thickness h. 

Thus, for the each TM mode of the hyperbolic waveguide two cut-off frequencies exist: b{Vc‘^m) = bo and 6(V'i m) = 0. 
For each TM mode of conventional dielectric waveguide only single cut-off frequency exists Vc^m. 
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C. Poynting vector of the guided waves in the hyperbolic waveguide 

The Poynting vector defines density of the radiation energy flux and direction of wave’s energy propagation. It is 
instructive to consider an averaged projection of the Poynting vector along the OZ axis. For TE wave it can be found 
from the equation 




(16) 


and for TM wave it follows from the relation 


{s;) = —{EiHy + E,H;). 


(17) 


The relations (16) and 0 can be rewritten with taking into account equation © for the TE wave and equation © 
for the TM wave as 


{S.)tE = i-^\Ey\\ {S.)tM = 

oTT 1^1 j oTT Sj 


(18) 


where j = 1,i, 3 for fij and j = for Sj depending on the layer under consideration. 

As was obtained in the previous two subsections, the effective index of refraction Uef for TE and TM guided modes 
in the hyperbolic waveguide can achieve null value. In these cases the averaged energy flux along the guided wave 
propagation axis OZ will be zero. Thus, the effect of slowing light in the waveguide takes place in these cases. 


IV. CONCLUSION 

The special kind of the hyperbolic slab waveguide is considered here. In the case of < 0 and 5e > 0 the modes of 
directed waves were found. The TE wave is ordinary wave, but TM wave is extraordinary in the waveguide under 
consideration. If > 0 and Se < 0 TM wave is not confined, and TE guided waves are identical with waves in a 
conventional waveguide. It is the reason to study the case of hyperbolic waveguide with < 0 and Se > 0. 

The effective index of refraction for TE wave obeys following inequality 0 < Uef < where rii is the index of 
refraction of the waveguide’s core. In the case of TM wave the effective index of refraction varies within the limits 
0 < rief < Tie, where Ue is the index of refraction of the extraordinary wave in hyperbolic medium. (It is assumed that 
He < rii.) Thus, in both cases the effective index of refraction can be equal to zero, that leads to great phase velocity 
of the wave and zero Poynting vector component along the wave’s propagation axis in the hyperbolic waveguide 
discussed above. Hence the light wave can be slowed down in this waveguide. In the case of conventional (elliptic) 
anisotropic dielectric waveguide the effective index of refraction lies in the range Uq < Uef < rii foi* TM wave, and in 
the range tIq < ri^f < rii for TE wave. 

The dispersion relations for the case of TE and TM waves are derived. It was shown that for the TM wave the 
number of guided modes is limited. Each of these modes have two cut-off frequencies. One of them corresponds to 
mode appearance, another corresponds to mode disappearance. There is region of parameters in which the only single 
mode exists in this waveguide. It is worth noting that this phenomenon is unavailable in the case of a conventional 
waveguide. Usually the number of modes increases with core thickness, and only single cut-off frequency exists. 
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